Determinantal point processes 
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Abstract 

We present a list of algebraic, combinatorial, and analytic mechanisms 
that give rise to determinantal point processes. 



1 Introduction 

Let X be a discrete space. A (simple) random point process V on X is a 
probability measure on the set 2 X of all subsets of X. V is called determinantal 
if there exists a \X\ x \X\ matrix K with rows and columns marked by elements 
of X, such that for any finite Y = (yi, . . . , y n ) C X one has 



Pr{X G 2 X | Y C X} = det[K{ yi , Vj ) 



A similar definition can be given for X being any reasonable space; then the 
measure lives on locally finite subsets of X. 

Determinantal point processes (with X = R) have been used in random 
matrix theory since early 60's. As a separate class determinantal processes 
were first singled out in |Mac75| to model fermions in thermal equilibrium, cf. 
[Ben73j . and the term 'fermion' point processes was used. The term 'deter- 
minantal' was introduced in [BorOOa], for the reason that the particles of the 
process studied there were of two kinds; particles of the same kind repelled, 
while particles of different kinds attracted. Nowadays, the expression 'determi- 
nantal point process (or field)' is standard. 

There are several excellent surveys of the subject available, see [SosOOj, 
|Lyo03| , |Joh05| , |Kon05| , [Hou06j , [Sos06j . The reader may find there a detailed 
discussion of probabilistic properties of determinantal processes as well as a 
wide array of their applications; many applications are also described in various 
chapters of this volume. 
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The goal of the present note is to bring together all known algebraic, combi- 
natorial, and analytic mechanisms that produce determinantal processes. Many 
of the well-known determinantal processes fit into more than one class described 
below. However, none of the classes is superseded by any other. 

2 Generalities 

Let X be a locally compact separable topological space. A point configuration 
X in X is a locally finite collection of points of the space X. Any such point con- 
figuration is either finite or infinite. For our purposes it suffices to assume that 
the points of X are always pairwise distinct. The set of all point configurations 
in X will be denoted as Conf(X). 

A relatively compact Borel subset A C X is called a window. For a window A 
and X E Conf(X), set Na(X) = \ AnX\ (number of points of X in the window). 
Thus, Na can be viewed as a function on Conf(X). We equip Conf(X) with the 
Borel structure generated by functions Na for all windows A. 

A random point process on X is a probability measure on Conf(X). 

Given a random point process, one can usually define a sequence {p n }^ =1 , 
where p n is a symmetric measure on X n called the nth correlation measure. 
Under mild conditions on the point process, the correlation measures exist and 
determine the process uniquely, cf. |Len73j . 

The correlation measures are characterized by the following property: For 
any n > 1 and a compactly supported bounded Borel function / on X™ one has 



where the sum on the right is taken over all n-tuples of pairwise distinct points 
of the random point configuration X. 

Often one has a natural measure ^onl (called the reference measure) such 
that the correlation measures have densities with respect to p® n , n = 1,2,... . 
Then the density of p n is called the nth correlation function and it is usually 
denoted by the same symbol "pn" • 

If X C R and p is absolutely continuous with respect to the Lebesgue 
measure, then the probabilistic meaning of the nth correlation function is that 
of the density of probability to find an eigenvalue in each of the infinitesimal 
intervals around points x%, x%, . . . x n : 

p n (xi,x 2 , ■ ..x n )p{dxi) ■ ■■ p(dx n ) 




(2.1) 



Pr {there is a particle in each interval (xi,Xi + dxi)}. 
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On the other hand, if fx is supported by a discrete set of points, then 

p n (x 1 ,x 2 , ■ ■ ■ x n )n(xi) ■ ■ ■ fj,(x n ) 

= Prjthere is a particle at each of the points Xi}. 

Assume that we are given a point process V and a reference measure such 
that all correlation functions exist. The process V is called determinantal if 
there exists a function J(:IxI-tC such that 

p n (x 1 , . . . ,x n ) = det[K(xi,Xj)]? j=1 , n = l,2, .... (2.2) 

The function K is called a correlation kernel of V . 

The determinantal form of the correlation functions (|2.2p implies that many 
natural observables for V can be expressed via the kernel K. We mention a few 
of them. For the sake of simplicity, we assume that the state space X is discrete 
and fi is the counting measure; under appropriate assumptions, the statements 
are easily carried over to more general state spaces. 

• Let I be a (possibly infinite) subset of X. Denote by Kj the operator 
in £ 2 (I) obtained by restricting the kernel K to /. Assume that Kj is a 
trace class operator. Then the intersection of the random configuration 
X with / is finite almost surely and 

Pr{| X nI\=N} = ^ det (l - zKj) ^ . 

In particular, the probability that X n / is empty is equal to 
Pr{X n / = 0} = det (l - Kj 



More generally, if Ji, . . . , I m is a finite family of pairwise nonintersecting 
intervals such that the operators Kj 1 , . . . , Kj m are trace class then 

Pr{|An/i| = A^i,...,|Xn/ m | = N m } 



ffiiW dz?K..dz. 



det ( 1 - z 1 K h z m K Ir , 



III 



-Z m = l 

(2.3) 



• Slightly more generally, let ^ be a function on X such that the kernel 
(1 — 4>(x))K(x,y) defines a trace class operator (1 — 4>)K in f 2 (X). Then 




E [I <P(xi) = det(l - (1 - 0)1T). (2.4) 
Specifying 4> = Sj=i(l — z j)^-ij l ea -ds to (|2.3p . 



1 For discrete X, a convenient sufficient condition for Ki to be of trace class is 
T, x , y£ i \K{x,y)\ < oo. 
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• For Jcl such that Kj is trace class and det(l — Kj) ^ 0, and arbitrary 
pairwise distinct locations {x±, . . . , x n } C /, n = 1, 2, . . . , set 

Ji. 

n (x\, . . . , x n ) — Pr {there is a particle at each of the points X{ 

and there are no other particles in /}. 

These are sometimes called Janossy measures. One has 

Ji,n(xi,...,Xn) = det(l - Ki) ■ det[L I (x i ,x j )]l j=1 , (2.5) 

where Lj is the matrix of the operator Kj(l — Kj)^ 1 . 

Simple linear-algebraic proofs of (|2.4p and (|2.5p can extracted from the proof 
of Proposition A. 6 in [BorOOb] . We also refer to Chapter 4 in this volume for a 
detailed discussion of (12.3p - (12.5p and many related identities. 

3 Loop-free Markov chains 

Let X be a discrete space, and let P = [P X y]x,yeX be the matrix of transition 
probabilities for a discrete time Markov chain on X. That is, P xy > for all 
x, y S X and 



Let us assume that our Markov chain is loop-free, i.e. the trajectories of the 
Markov chain do not pass through the same point twice almost surely. In other 
words, we assume that 



This condition guarantees the finiteness of the matrix elements of the matrix 



Indeed, (P k ) X y is the probability that the trajectory started at x is at y after 
kth step. Hence, Q xy is the probability that the trajectory started at x passes 
through y ^ x, and since there are no loops we have Q xy < 1. Clearly, Q xx = 0. 
The following (simple) fact was proved in [Bor08aJ. 

Theorem 3.1 For any probability measure ir = [tt x ] X £X on consider the 
Markov chain with initial distribution ir and transition matrix P as a probability 
measure on trajectories viewed as subsets of X. Then this measure on 2^ is a 
determinantal point process on X with correlation kernel 




P xy = 1 for any x £ X. 



(P k ) xx = for any k > and x <E X. 



Q = P + P 2 + P 3 + ... . 



K(x,y) 



k x + {ttQ) x - Q 



yx- 
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Note that the correlation kernel is usually not self-adjointQ and self-adjoint 
examples should be viewed as "exotic" . One such example goes back to |Mac75j . 
see also §2.4 of [SosOOj : It is a 2-parameter family of renewal processes — pro- 
cesses on Z or M with positive i.i.d. increments. Theorem 13 . 1 1 implies that if we 
do not insist on self-adjointness then any process with positive i.i.d. increments 
is determinantal. 

4 Measures given by products of determinants 

Let X be a finite set and N be any natural number no greater than \X\. Let <3? ra 
and Vl/ n , n = 1, 2, . . . , N, be arbitrary complex-valued functions on X. To any 
point configuration X £ Conf (X) we assign its weight W{X) as follows: If the 
number of points in X is not ./V then W(X) = 0. Otherwise, using the notation 
X = {xi, . . . , xn}, we have 

W(X) = det [Mxj))Z=i det [*i(*i)G=i • 
Assume that the partition function of our weights does not vanish 

XGConf(X) 

Then the normalized weights W(X) = W(X)/Z define a (generally speaking, 
complex valued) measure on Conf(X) of total mass 1. Such measures are called 
biorthogonal ensembles^ For complex valued point processes we use (12. 1|) to 
define their correlation functions. 

An especially important subclass of biorthogonal ensembles consists of or- 
thogonal polynomial ensembles, for which X must be a subset of C, and 

N 

w(x)= n \ Xl ~ Xj \ 2 - Hw(xi) 

l<i<j<N i=l 

for a function w : X — > M+, see e.g. [K6n05] and Chapter 4 of this volume. 

Theorem 4.1 Any biorthogonal ensemble is a determinantal point process. Its 
correlation kernel has the form 

N 

K{x,y)=Y J [G-%*iW*j(v), 

where G = [Gyf^- =1 is the Gram matrix: G^ = YlxeX ®i(. x )^j( x )^ 

2 In fact, it can be written as a sum of a nilpotent matrix and a matrix of rank 1. 

3 This term was introduced in |Bor99] and is now widely used. 

4 The invertibility of the Gram matrix is implied by the assumption Z 7^ 0. 
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The statement immediately carries over to X being an arbitrary state space 
with reference measure /x; then one has Gij = j x ^i(x)^j(x)n(dx). 

Probably the first appearance of Theorem 14.11 is in the seminal work of 
F. J. Dyson [Dys62a|, where it was used to evaluate the correlation functions 
of the eigenvalues of the Haar-distributed N x N unitary matrix. In that case, 
X is the unit circle, \x is the Lebesgue measure on it, 



5«-l 



1, 



1. 



,N, 



and the Gram matrix G coincides with the identity matrix. 

In the same volume, Dyson |Dys62b| introduced a Brownian motion model 
for the eigenvalues of random matrices (currently known as the Dyson Brownian 
motion), and it took more than three decades to find a determinantal formula 
for the time-dependent correlations of eigenvalues in the unitarily invariant case. 
The corresponding claim has a variety of applications; let us state it. Again, 
for simplicity of notation, we work with finite state spaces. 

Let . . . , £( fc ) be finite sets. Set X = £W U • • • U £(*). Fix a natural 
number N. Let 



$i : X^ -f C, : X^ 

T jJ+1 : X (j) x X (j+1) -» C, 



1,. 

,fc 



N 



be arbitrary functions. To any X G Conf(X) assign its weight W(X) as follows. 
If X has exactly iV points in each X^\ j = 1, . . . , k then denoting X n 3£V' = 
. . . , xjy } we have 



WpQ = det 



v 



det 



A' 



x det 



T k -\,k{x 



v 



det 



v 



(4.1) 



otherwise W(X) = 0. 

As for biorthogonal ensembles above, we assume that the partition func- 
tion of these weights is nonzero and define the corresponding normalized set 
of weights. This gives a (generally speaking, complex valued) random point 
process on X. 

In what follows we use the notation 

Z X 

(hi * f){y) = hx(x)f(x, y), (g * h 2 )(x) = ^ g(x, y)h 2 (y) 



for arbitrary functions f(x,y), g(x,y), hi(x), h,2{x), where the sums are taken 
over all possible values of the summation variables. 
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Theorem 4.2 The random point process defined by J^, l\j is determinantal. 
The correlation kernel on x 3£( q >, p,q = 1,...,N, can be written in the 
form 

K{x<r\yM) = -l p>g • (T M+1 * • • • * V 1)P )(^,xW) 

N 

+ E l G ~% *n,2*-~* W (^ (p) ) CWi * • • • * r *-i,* * 

(4.2) 

where the Gram matrix G = [Gij] i - =1 is defined by 

G tj = * 71,2 * • • • * T fc _i ifc * i, j = 1, . . . , N. 

Similarly to Theorem 14.11 the statement is easily carried over to general 
state spaces 

Theorem l4.2l is often referred to as the Eynard-Mehta theorem, it was proved 
in |Eyn98| and also independently in |Nag98| . Other proofs can be found in 
|Joh03j . |Tra04j . |Bor05j . 

The algebraically "nice" case of the Eynard-Mehta theorem, which e.g. takes 
place for the Dyson Brownian motion, consists in the existence of an orthonor- 
mal basis in each L 2 (X^), j = 1, . . . , k, such that 



T j>j+1 (x, y)=Y, c i,i+i;* j = 1, 2, . . . , k - 1, 



i>l 

for some constants Cjj+i-j, and 

Span{sW, . . . ,3$} = Span^j, . . . , $ N }, 
Span{E<* \ . . . , E%>} = Span{*i, . . . , ^^v}. 

Then, with the notation c fe)J;i = c k}k +i-,iCk+i,k+2;i • ■ ■ Q-l,i;i> flM]) reads 

( N 



K{x {p \y {q) ) 



i>7V 



The ubiquitousness of the Eynard-Mehta theorem in applications is ex- 
plained by the combinatorial statement known as the Lindstrom-Gessel-Viennot 
(LGV) theorem, see | Ste90] and references therein, that we now describe. 

Consider a finitqj directed acyclic graph and denote by V and E the sets 
of its vertices and edges. Let w : E — > C be an arbitrary weight function. For 



3 The assumption of finiteness is not necessary as long as the sums in (|4.3f) converge. 
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any path it denote by w(ir) the product of weights over the edges in the path: 
w(ir) = Yl e£7T w(e). Define the weight of a collection of paths as the product of 
weights of the paths in the collection (we will use the same letter w to denote 
it). We say that two paths n\ and tt2 do not intersect (notation 7ri PI iT2 = 0) 
if they have no common vertices. 

For any u, v G V, let H(u, v) be the set of all (directed) paths from u to v. 

Set 

T(u,v)= ( 43 ) 

Theorem 4.3 Let (u\, . . . ,u n ) and (v\, . . . ,v n ) be two n-tuples of vertices of 
our graph, and assume that for any nonidentical permutation a E S(n), 

{(tT%, ... ,7T n ) | VTj G II (ui,V a (i)) , 7Tj D TVj = 0, i, j = 1, . . . , Tl} = 0. 

Then 

^2 w(7ri,...,7r n ) = det [T(ui,Vj)]^ j=1 . 

TTl&U(ui,Vl),...,TT n eU(u n ,V n ) 

iTiniTj=0, i,j=l,...,n 

Theorem 14.31 means that if, in a suitable weighted oriented graph, we have 
nonintersecting paths with fixed starting and ending vertices, then the distribu- 
tions of the intersection points of these paths with any chosen "sections" have 
the same structure as (|4.ip . and thus by Theorem l4.2l we obtain a determinantal 
point process. 

A continuous time analog of Theorem 14.31 goes back to [Kar59| , who in 
particular proved the following statement (the next paragraph is essentially a 
quotation). 

Consider a stationary stochastic process whose state space is an interval on 
the extended real line. Assume that the process has strong Markov property 
and that its paths are continuous everywhere. Take n points x\ < ■ ■ ■ < x n 
and n Borel sets E\ < ■ ■ ■ < E n , and suppose n labeled particles start at 
xi, . . . ,x n and execute the process simultaneously and independently. Then the 
determinant det [Pt(xi, Ej)]™- V with P t (x,E) being the transition probability 
of the process, is equal to the probability that at time t the particles will be 
found in sets E\,... ,E n respectively without any of them ever having been 
coincident in the intervening time. 

Similarly to Theorem 14.31 this statement coupled with Theorem 14.21 leads 
to determinantal processes, and this is exactly the approach that allows one 
to compute the time-dependent eigenvalue correlations of the Dyson Brownian 
motion. 

We conclude this section with a generalization of the Eynard-Mehta theorem 
that allows the number of particles to vary. 



S 



Let %x, . . . , 3Cn be finite sets, and 

0n( * > " ) ■ 3£n—l X 3C n 

n (virt, • ) : 3t n 



n 
n 



2,...,N, 
l,...,N, 



1, 



^•(•):3Ejv-»C, j 

be arbitrary functions on the corresponding sets. Here the symbol virt stands for 
a "virtual" variable, which is convenient to introduce for notational purposes. 
In applications, virt can sometimes be replaced by +00 or —00. 
Let c(l), . . . ,c(N) be arbitrary nonnegative integers, and let 



< t 



N 



,N-1 < 



< t 



N-l 



,N-2 < 



< tt 



*0< 



< t; 



^ ^ »c(N) — "0 - ' ' ' - fc c (Ar-l) - "0 - ' ' ' - b c{2) ~ "0 ^ " " " ^ u c{l) 

be real numbers. In applications, these numbers may refer to time moments of 
an associated Markov process. Finally, let 

' ' ' ) : % n x ~* C ' n = l,...,iV, a = 1, . . . ,c(n), 

be arbitrary functions. 

Set X = (Xi U ■ ■ ■ U Xj) U • • • U (3Ln U • • • U 3Ln) with c(n) + 1 copies of each 
3tJ§, and to any X € Conf(X) assign its weight W(X) as follows. 

The weight is zero unless X has exactly n points in each copy of X n , 

n = 1, . . . , N. In the latter case, denote the points of X in the mth copy of X n 
by x^(tj^), k = 1, . . . , n, and set 



n=l 



det 

c(n) 



„n— 1 /j.n— 1 



(*c(n))) 



fcj=l 



a=l 



1 « 
fcJ=l 



det 



N 



k,l=l 

(4.4) 

where 1 ( ■ ) = virt for all n = 1, . . . ,N. 

Once again, we assume that the partition function does not vanish, and 
normalizing the weights we obtain a (generally speaking, complex valued) point 
process on X. 

We need more notation. For any n = 1,...,N and two time moments 



*™ > tb we define 



T+n 4-n = T+n 4-n # 1+n +n 



u b+l> b 



rj-n 



T-fn 4-n . 
c(i)' ° 



For any time moments t" 1 > with (oi,ni) 7^ (02,^2)5 we denote the convo- 
lution over all the transitions between them by (f)^ a i > a 2 > ; 



4> (ta i ' ta 2 ) = r t «i t ni * <^ ni+1 * r ni+1 * • • • * n2 * 7^. 



(«2) 



Instead of c(n) + 1 copies of £„ one can take same number of different spaces, and a 
similar result will hold. We decided not to do it in order not to clutter the notation anymore. 
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If there are no such transitions, i. e. if t™\ < tfj® or (ai,rii) = (a2,n2), we set 
Furthermore, define the "Gram matrix" G = [Gfc/k j =1 by 

and set 

$rj™ = </>(C*o0 i = 1, . . . , JV. 

Theorem 4.4 TTie random point process on X defined by ^4-4\ ) is determinan- 
tal. Its correlation kernel can be written in the form 

K(Q, Xl ;Q,x 2 ) = -cj>^l\x 2 ,x x ) 

711 N "2 

+ EE * ^'^^(virt^i) vl/*" 2 (x 2 ). 

i=i j=i 

One proof of Theorem l4,4l was given in [Bor08b| : another proof can be found 
in Section 4.4 of |F or08| . Although we stated Theorem 14.41 for the case when 
all sets X n are finite, one easily extends it to a more general setting. 



5 L-ensembles 

The definition of L-ensembles is closely related to (|2.5p . 

Let X be a finite set. Let L be a |X| x |X| matrix whose rows and column 
are parameterized by points of X. For any subset X C X we will denote by Lx 
the symmetric submatrix of L corresponding to X: Lx = [L(xi,Xj)] x , x . eX ■ 
If determinants of all such submatrices are nonnegative (e.g., if L is positive 
definite), one can define a random point process on X by 

1 1 det(l + L) 

This process is called the L-ensemble. 

The following statement goes back to [Mac75] . 

Theorem 5.1 The L-ensemble as defined above is a determinantal point pro- 
cess with the correlation kernel K given by K = L(l + L) -1 . 

Take a nonempty subset 2} of X and, given an L-ensemble on X, define a 
new random point process on 2) by considering the intersections of the random 
point configurations X C X of the L-ensemble with 2J> provided that these 



10 



point configurations contain the complement 2) of 2) in X. It is not hard to see 
that this new process can be defined by 



Pr{T} 




Y G Conf(2J). 



det(l«g + L) 



Here lm is the block matrix 



1 





where the blocks correspond to the splitting 



X = 2) U 2). This new process is called the conditional L-ensemble. The next 
statement was proved in [Bor05]. 

Theorem 5.2 The conditional L-ensemble is a determinantal point process 
with the correlation kernel given by 



Note that for 2) = X, Theorem 15,21 coincides with Theorem 15.11 

Not every determinantal process is an L-ensemble; for example, the pro- 
cesses afforded by Theorem 14.11 have exactly N particles, which is not possible 
for an L-ensemble. However, as shown in [Bor05j, every determinantal process 
(on a finite set) is a conditional L-ensemble. 

The definition of L-ensembles and Theorems 15.11 15.21 can be carried over 
to infinite state spaces X, given that L satisfies appropriate conditions. In 
particular, the Fredholm determinant det(lw + L) needs to be well defined. 

Although L-ensembles do arise naturally, see e. g. [BorOOb] . they also con- 
stitute a convenient computation tool. For example, proofs of Theorems 14.21 
and 14.41 given in [Bor05j and [Bor08b] represent the processes in question as 
conditional L-ensembles and employ Theorem 15. 2[ 

Here is another application of Theorem 15.21 

A random point process on (a segment of) Z is called one-dependent if for 
any two finite sets A, B C Z with dist(A, B) > 2, the correlation function 
factorizes: p\ A \+\B\(A U B) = p\ A \ (A)p\ B \ (B). 

Theorem 5.3 Any one-dependent point process on (a segment of) Z is deter- 
minantal. Its correlation kernel can be written in the form 



K = la, - (la, + L) 



i 



x — y > 2 
x — y = 1 



K(x,y) 



= < 



-1, 



y—x+Y 



E (-i) 



r-l 



E 



R l ,h R h,h • ' ' R lr-l,lr X <V, 



r=l 



x=lo<li<-<l r =y+l 



where R a , h = p b - a (a,a + 1, ... ,6 - 1). 



Details and applications can be found in [Bor09]. 
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6 Fock space 



A general construction of determinantal point processes via the Fock space 
formalism can be quite technical, see e. g. [Lyt02| , so we will consider a much 
simpler (however nontrivial) example instead. 

Recall that a partition A = (Ai > A2 > • • • > 0) is a weakly decreasing 
sequence of nonnegative integers with finitely many nonzero terms. We will use 
standard notations |A| = Ai + A2 + • • • for the size of the partition and £(\) for 
the number of its nonzero parts. 

The poissonized Plancherel measure on partitions is defined by 

Pr{A} = e- 2 ( ni^^- '- A j +j) glx Y 

where 8 > is a parameter, and L is an arbitrary integer > £(X). We refer to 
[BorOOb] . [JohOl] and references therein for details. 

It is convenient to parameterize partitions by subsets of Z' := Z + |: 

A^£(A) = {Aj-i + ^cZ'. 

The pushforward of (|6.1[) via C defines a point process on Z', and we aim to 
show that it is determinantal. We follow [OkoOl]; other proofs of this fact can 
be found in [BorOOb] . [JohOlj . 

Let V be a linear space with basis {A; | k £ Z'}. The linear space A^V is, 
by definition, spanned by vectors 

^ = £1 AS2AS3A... , 

where S = {s± > S2 > ■ ■ ■ } C Z' is such that both sets S + = S \ Z' <0 and 
S- = Z' <0 \ S are finite. We equip \~2~V with the inner product in which the 
basis {vs} is orthonormal. 

00 

Creation and annihilation operators in A~z~V are introduced as follows. The 
creation operator tp^ is the exterior multiplication by k: tp^ (/) = k A / . The 
annihilation operator is its adjoint. These operators satisfy the canonical 
anti-commutation relations 

i>ki>i + i>*i>k = h,l , k,l G Z'. 

Observe that 

jv s , keS, 

^ 5 = \o, k{s. (6 - 2) 

Let C be the charge operator: Cvs = (|£+| — |£_|)us. One easily sees that 
the zero-charge subspace kerC C A~V is spanned by the vectors vc(\) with A 
varying over all partitions. The vacuum vector 

^Vac — 2^ 2^ 2 ^ ' ' ' 
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corresponds to the partition with no nonzero parts. 

Define the operators a n = ^fcez' ^k-n^ti n G Z \ {0}. Although the sums 
are infinite, the application of a n to any vs yields a finite linear combination 
of basis vectors. These operators satisfy the Heisenberg commutation relations 

a m a n - a n a m = m 8 n - m , m,n eZ\{0}. 

For any 8 > 0, define T±(0) = exp(#a±i). It is not difficult to show that 

T* ± (6) = T T (6), r+(0)r_(0') = e ee ' ■ V-(9')r + (6), T + (9)v vac = v vac . 

(6.3) 

One also proves that 

r_(e)v vac - 2^ L — — — v l v c{X) , 

A V n 4 =i(Ai -i + L)\ J 

where the sum is taken over all partitions, cf. (|6,ip . This implies, together with 
(|6.2p . that for any n > 1 and xi, ■ ■ ■ ,x n € Z', the correlation function of our 
point process can be written matrix element 



p n (xi,... ,x n ) = e 02 ( (II ^xiipxij r -(#) Wvac,r_(6»)t; vac J . 
Using ()6.3p we obtain 

xa ^Vac j ^vac (6-4) 



\i=l 

where 



^ k = G^ k G-\ n = or k G~\ G = r + {e)r_(9)- 1 



Theorem 6.1 We have 

p n (xi,...,x n ) =det[K(xi,Xj)]" j=1 , (6.5) 

where K(x,y) = (^ x ^l v vac , v V3C ) . 

The passage from (|6.4p to (|6.5p is an instance of the fermionic Wick theo- 
rem; it uses the fact that ^f x and \&* are linear combinations of ^'s and ^*'s 
respectively, together with the canonical anti-commutation relations. 

A further computation gives an explicit formula for the correlation kernel: 

K{x,y) = 2 y+2 _ +2 y 2 = ]T J x+k J y+k , 
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where Jk = </fc(20) are the J-Bessel functions. This is the so-called discrete 
Bessel kernel that was first obtained in [BorOOb], [JohOlj . 

We refer to [Qko01| and [Oko03| for far-reaching generalizations of Theorem 
16.11 and to |Lyt02| for a general construction of determinantal processes via 
representations of the canonical anti-commutation relations corresponding to 
the quasi-free states. 

7 Dimer models 

Consider a finite planar graph Q. Let us assume that the graph is bipartite, i. e. 
its vertices can be colored black and white so that each edge connects vertices 
of different colors. Let us fix such a coloring and denote by B and W the sets 
of black and white vertices. 

A dimer covering or a domino tiling or a perfect matching of a graph is a 
subset of edges that covers every vertex exactly once. Clearly, in order for the 
set of dimer coverings of Q to be nonempty, we must have \B\ = \W\. 

A Kasteleyn weighting of Q is a choice of sign for each edge with the property 
that each face with mod 4 edges has an odd number of minus signs, and each 
face with 2 mod 4 edges has an even number of minus signs. It is not hard to 
show that a Kasteleyn weighting of Q always exists (here it is essential that the 
graph is planar), and that any two Kasteleyn weightings can be obtained one 
from the other by a sequence of multiplications of all edges at a vertex by —1. 

A Kasteleyn matrix of Q is a signed adjacency matrix of Q. More exactly, 
given a Kasteleyn weighting of Q, define a \B\ x \ W\ matrix & with rows marked 
by elements of B and columns marked by elements of W, by setting R(b, w) = 
if b and w are not joined by an edge, and R(b, w) = ±1 otherwise, where ± is 
chosen according to the weighting. 

It is a result of [Tem61| . [Kas67] that the number of dimer coverings of Q 
equals | det^|. Thus, if there is at least one perfect matching, the matrix R is 
invertible. 

Assume that detfi. ^ 0. Define a matrix K with rows and columns parame- 
terized by the edges of Q as follows: K(e, e') = &~ l (w, b'), where w is the white 
vertex on the edge e, and b' is the black vertex on the edge e' . The next claim 
follows from the results of |Ken97] . 

Theorem 7.1 Consider the uniform measure on the dimer covers of Q as a 
random point process on the set X of edges of Q. Then this process is determi- 
nantal, and its correlation kernel is the matrix K introduced above. 

The theory of random dimer covers is a deep and beautiful subject that has 
been actively developing over the last 15 years. We refer the reader to |Ken08] 
and references therein for further details. 
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8 Uniform spanning trees 



Let G be a finite connected graph. A spanning tree of Q is a subset of edges of 
Q that has no loops, and such that every two vertices of Q can be connected 
within this subset. 

Clearly, the set of spanning trees is nonempty and finite. Any probability 
measure on this set can be viewed as a random point process on the set X of 
edges of Q. We are interested in the uniform measure on the set of spanning 
trees, and we denote the corresponding process on X by V. 

Let us fix an orientation of all the edges of Q. For any two edges e = xy 
and / denote by K (e, /) the expected number of passages through /, counted 
with a sign, of a random walk started at x and stopped when it hits y. 

The quantity K (e, /) also has an interpretation in terms of electric networks. 
Consider Q with the fixed orientation of the edges as an electric network with 
each edge having unit conductance. Then K (e, /) is the amount of current 
flowing through the edge / when a battery is connected to the endpoints x and 
y of e, and the voltage is such that unit current is flowing from y to x. For this 
reason, K is called the transfer current matrix. 

This matrix also has a linear-algebraic definition. To any vertex v of Q we 
associate a vector a(v) € i 2 {X) (recall that X is the set of edges) as follows: 

1, if v is the tail of e, 
— 1, if v is the head of e, 
0, otherwise. 

Then [K (e, f)] e ,feX is the matrix of the orthogonal projection operator with 
image Span{a(i>)}, where v varies over all vertices of G, see c. g. [BcnOl . 

Theorem 8.1 V is a determinantal point process with correlation kernel K . 

Theorem 18. II was proved in |Bur93| : another proof can be found in [Ben01| . 
The formulas for the first and second correlation functions via K go back to 
[Kir 1847] and |Bro40] . respectively. We refer to |Lyo03| and references therein 
for further developments of the subject. 

Note that for planar graphs, the study of the uniform spanning trees may be 
reduced to that of dimer models on related graphs and vice versa, see |Bur93j . 
[KenOOj . |Ken04j . 

9 Hermit ian correlation kernels 

Let X be M. d or 7L d with the Lebesgue or the counting measure as the reference 
measure. Let K be a nonnegative operator in L 2 (xjl|. In the case X = R d , 

7 An Hermitian K must be nonnegative as we want det[K(xi,Xj)] > 0. 



eeX 
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we also require K to be locally trace class, i. e. for any compact B C 3t, the 
operator K ■ 1b is trace class. Then Lemma 2 in [SosOO shows that one can 
choose an integral kernel K(x,y) of K so that 

Trace(K-lB) fc = / K(xi, x 2 )K (aj 2 , x s ) ■ ■ ■ K(xk, xi)dxi ■ ■ ■ dxk, fc = l,2, ... 
For X = Z d , K(x, y) is just the matrix of K. 

Theorem 9.1 There exists a determinantal point process on X with the cor- 
relation kernel K(x,y) if and only if < K < 1, i. e. both K and 1 — K are 
nonnegative. 

Theorem 19.11 was proved in [SosOO] ; an incomplete argument was also given 
in |Mac75j . Remarkably, it remains the only known characterization of a broad 
class of kernels that yield determinantal point processes. 

Although only Theorem 14 . 1 1 with $j = ^j, Theorem 16. 1( and Theorem 18.11 
from the previous sections yield manifestly nonnegative kernels, determinantal 
processes with such kernels are extremely important, and they are also the 
easiest to analyze asymptotically, cf. [Hou06j. 

Let us write down the correlation kernels for the two most widely known 
determinantal point processes; they both fall into the class afforded by Theorem 

ECU 

The sine process on R corresponds to the sine kernel 

K siaC (x,y) = Sin < X ~y) = f~ 2 e 2inr Xe -2inry dT _ 

n{x - y) J_i 

The Fourier transform of the corresponding integral operator K sme in L 2 (R) 
is the operator of multiplication by an indicator function of an interval; hence 
K sme is a self-adjoint projection operator. 

The Airy point process^ on R is defined by the Airy kernel 

r^Airv, \ Ai(x)Ai'(y) -Ai'(x)Ai(y) f +oc ... 

K AlTy (x, y) = — — — — — — = / Ai(x + r)Ai(y + r)dr, 

x-y Jo 

where Ai{x) stands for the classical Airy function. The integral operator K AlTy 
can be viewed as a spectral projection operator for the differential operator 
■j-j — x that has the shifted Airy functions {Ai(x + T)} Tg R as the (generalized) 
eigenf unctions. 



Not to be confused with the Airy process that describes the time evolution of the top 
particle of the Airy point process, see Chapter 37 of the present volume. 
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10 Pfaffian point processes 



A random point process on X is called Pfaffian if there exists a 2 x 2 matrix 
valued skew-symmetric kernel K on X such that the correlation functions of 
the process have the form 

p n (xi,...,x n ) = ¥i[K(xi,Xj)}™ j=l , xi,...,x n e X, n = l,2, ... 

The notation Pf in the right-hand side stands for the Pfaffian, and we refer to 
[deB55] for a concise introduction to Pffafians. 

Pfaffian processes are significantly harder to study than determinantal ones. 
Let us list some Pffafian analogs of the statements from the previous sections. 

• A Pfaffian analog of the Fredholm determinant formula (|2.4p for the gen- 
erating functional can be found in Section 8 of RaiOOj. 

• A Pfaffian analog of the Eynard-Mehta theorem is available in [Bor05j. 

• Pfaffians can be used to enumerate nonintersecting paths with free end- 
points, sec [Ste90j. This leads to combinatorial examples for the Pfaffian 
Eynard-Mehta theorem. 

• Pfaffian L-ensembles and conditional L-ensembles are treated in [Bor05j. 

• Fermionic Fock space computations leading to a Pfaffian point process 
were performed in [Fer04| and |Vul07j . 

• Pfaffians arise in the enumeration of dimer covers of planar graphs that 
are not necessarily bipartite, see [Tem61| . [Kas67| . 
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